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Abstract. We present very accurate individual masses 
(1.2% relative accuracy) for the two components of 
Gl 570BC, an interferometric and double-lined spectro- 
scopic binary system. They were obtained from new high 
accuracy radial velocity and angular separation measure- 
ments, analysed together with previously published mea- 
surements. From those data we determine a much im- 
proved orbit through a simultaneous least square fit to 
the radial velocity, visual, and parallax information. The 
derived masses and absolute magnitudes generally vali- 
date the theoretical and empirical mass-luminosity rela- 
tions around 0.5 Mg, but point towards some low level 
discrepancies at the 0.1 to 0.15 magnitude level. Forth- 
coming results of this observing program will extend the 
comparison to much lower masses with similar accuracy. 

Key words: stars: binaries - stars: low mass, brown 
dwarfs - stars: individual: Gl 570 - techniques: radial ve- 
locity - techniques: adaptive optics 



1. Introduction 

Fundamental determination of stellar masses from binary 
orbits is a most classical astrophysical discipline, last com- 
prehensively reviewed by Andersen (1991; 1998). Besides 
defining a mean mass-luminosity relation, used at many 
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places in astronomy to approximately convert the observ- 
able stellar light to the underlying mass (for instance to 
derive an initial mass function), accurate stellar masses 
in multiple systems provide what is perhaps the most de- 
manding and fundamental test of stellar evolution theory 
(e.g. Andersen 1991). Mass, the basic input of evolutionary 
models, is directly measured and the models must repro- 
duce the effective temperatures and radii (or luminosities) 
of both components, for a single age and a single chemical 
composition. Given the strong dependence of all stellar 
parameters on mass, this discriminating diagnostic how- 
ever only shows its power for relative mass errors <l-2%. 
In practice this has up to now restricted its use to double- 
lined detached eclipsing binaries. These systems are how- 
ever relatively rare: only 44 have yielded masses accurate 
enough to be included in Andersen's (1991) critical com- 
pilation, and few have appeared in the literature since 
then. Tidally induced rotational mixing may in addition 
affect the evolution of the short period eclipsing systems, 
perhaps sufficiently that they are not completely represen- 
tative of isolated stars. More seriously however, detached 
eclipsing systems poorly fill some interesting areas in the 
HR diagram. 

The lower main sequence is one major region with 
very few known eclipsing systems, as a result of the 
strong observational and intrinsic biases against observing 
eclipses in faint and physically small stars. Only three well 
detached eclipsing binaries are known with significantly 
subsolar component masses: YY Gem (MOVe, 0.6+0.6 
M Q ; Bopp 1974; Leung & Schneider 1978), CM Dra 
(M4Ve, 0.2+0.2 M ; Lacy 1977, Metcalfe et al. 1996), and 
the recently identified GJ 2069A (M3.5Ve, 0.4+0.4 M ; 
Delfosse et al. 1999a). Detailed observational checks of 
evolutionary models (e.g. Paczynski & Sienkiewicz 1984; 
Chabrier & Baraffe 1995) have therefore heavily relied 
on the first two of these systems, even though they are 
in some respects non-ideal for this purpose: both bina- 
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ries have two nearly equal mass components, so that the 
strength of the differential comparison of two stars with 
different masses but otherwise equal parameters is largely 
lost. Also, all three are chromospherically very active, due 
to tidal synchronisation of their rotation with the short or- 
bital period. As a consequence, they may have untypical 
colours for their bolometric luminosity. 

Angularly resolved spectroscopic binaries provide stel- 
lar masses in parts of the HR diagram where eclipsing 
systems are missing, though to date these systems have 
not matched the ~1% accuracy which can be obtained in 
detached eclipsing systems. For M dwarfs in particular, 
the best representation to date of the empirical M-L rela- 
tion (Henry & Mc Carthy 1993, hereafter HMcC; Henry, 
Franz, Wasserman et al., 1999) still has to rely in part on 
some fairly noisy mass determinations. For several years 
(Perrier et al. 1992) we have therefore been following up 
with high angular resolution some low mass spectroscopic 
systems found with the CORAVEL or ELODIE radial 
velocity spectrographs. This follow-up initially used one- 
dimensional (ID) IR speckle, then two-dimensional (2D) 
IR speckle, and now uses adaptive optics imaging. As a 
progress report on this program and as an illustration of 
our methods, we present here much improved parameters 
for the double-lined interferometric binary Gl 570BC. The 
~1% accuracy for the 2 masses improves by an order of 
magnitude on our earlier measurements (Mariotti et al. 
1990) of this system, and is getting close to what is ob- 
tainable for eclipsing systems. 

The Gl 570 system comprises the V=5.7 K4V star 
Gl 570A (also HR 5568, HD 131977, HIP 73184, FK5 
1391), and at a projected distance of 25" the close 
Gl 570BC pair (also HD 131976, HIP 73182) which is 
the subject of the present paper. As discussed below, 
the orbital period of the close pair is fairly short, only 
10 months. Thanks to its small distance of only 6 pc it is 
nonetheless usually well resolved by the diffraction limit 
of 4m-class telescopes. The separation within Gl 570BC 
on the other hand always remains less than 0.2", so that 
all seeing-limited mesurements refer to integrated proper- 
ties of the close pair. Its integrated magnitude is V=8.09, 
and its joint spectral type is M1V (Henry et al. 1994; Reid 
et al. 1995). The three components have common paral- 
lax, proper motion, and radial velocity. They are therefore 
gravitationnally bound, with the projected separation of 
the wide pair (about 120 AU at the distance of the system) 
corresponding to P ~ 500 years. Formation of the system 
could either result from the fragmentation of a single gas 
clump or have involved some capture(s). The latter pro- 
cess however only remains efficient at the high densities 
which characterize very rich star forming clusters, whose 
lifetimes are very much shorter than the hydrogen burn- 
ing timescale in a K4 dwarf. For all practical purposes, the 
three stars can thus in both cases be considered as coeval 
and as formed from the same gas. 



2. Observations and results 

2.1. Angular separation measurements 

A variety of astrometric observations of different kinds and 
of different quality have been collected for Gliese 570BC 
over the years. The close pair was first resolved by Mari- 
otti et al. (1990) who obtained infrared ID speckle obser- 
vations on three occasions and derived a first visual orbit. 
We refer to their paper for the description of the instru- 
ment and the observing and data reduction procedures. 
HMcC later measured Gl 570BC once with a 2D infrared 
speckle instrument. For readers' convenience, we list these 
published measurements in Table |l| together with our new 
2D speckle and adaptive optics (hereafter AO) observa- 
tions. 

Two measurements were obtained in February 1991 
and April 1991 using the speckle mode of the 2D infrared 
imagers then installed respectively at the KPNO 4.2m 
telescope and the CFHT 3.6m telescope. Each imager was 
designed to permit acquisition of exposures short enough 
to substantially freeze the seeing under standard atmo- 
spheric conditions (t exp w 50 to 100 ms) in bands H and 
K, independently of the overhead due to read-out time 
and data transfer time. Several sequences of a few hundred 
such short exposures were obtained, alternating every few 
minutes between the source and a nearby unresolved star 
(usually Gl 570A) used as a point spread function (PSF) 
calibrator. The whole observation took one hour or less. 

In principle, this observing procedure allows an almost 
simultaneous PSF calibration, and consequently estimates 
the visibility modulus with 1 to 5% accuracy. We used a 
software package written specifically for this type of data 
reduction by E. Tessier (Tessier et al. 1994). It produces 
an unbiased visibility estimator for the source and then ex- 
tracts the binary parameters and their estimated variance 
from this visibility. The actual detector scale and position 
angle (P. A.) origin for those observations were calibrated 
from observations of the astrometric binary £ Aqr (Heintz 
1989). 

Most of the new measurements were obtained at 
the 3.6m Canada-France-Hawaii Telescope (CFHT) on 
top Mauna Kea, using the CFHT Adaptive Optics 
Bonnette (Arsenault et al. 1994, Rigaut et al. 1998) 
and two different infrared cameras (Nadeau et al. 
1994, Doyon et al. 1998). Delfosse et al. (1999b) pro- 
vide a detailed description of the observing procedure, 
which we therefore don't repeat here. Two observations 
were also obtained with the ESO 3.6m telescope (La Silla, 
Chile) AO system COME-ON+ (Rousset & Beuzit 1999), 
then ADONIS (Beuzit et al. 1999), equipped with the 
SHARP-II infrared camera, using a similar observing pro- 
cedure. This COME-ON+ observation has been published 
in Mariotti et al (1991). For recent CFHT measurements 
the corrected point spread function obtained from the AO 
system was synthesized from simultaneous recordings of 
the wavefront sensor measurements and deformable mir- 
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Table 1. Angular separation measurements 



Date 


Pos. Angle 


Separation 


Ref. 




degree 


arc second 




46961 





-0.060±0.010 


1 


46961 


90 


-0.090±0.010 


1 




4^ 

^±o 


-U.11U3ZU.U1U 


i 
i 


46961 


135 


-0.020±0.040 


1 


47071 





-0.154±0.025 


1 


47071 


90 


-0.097±0.025 


1 


47274 





-0.075±0.020 


1 


47274 


45 


-0.114±0.010 


1 


47274 


90 


-0.094±0.020 


1 


47343 





-0.130±0.010 


1 


47343 


45 


-0.173±0.010 


1 


47343 


90 


-0.114±0.020 


1 


Date 


p (arcsec) 


9 (degree) 


Ref. 


47934 


0.153±0.006 


220.0±2.0 


2 


48312 


0.194±0.010 


210.2±3.0 


3 


48372 


0.146±0.020 


203.0±5.0 


4 


49470 


0.159±0.004 


222.4±1.0 


5 


49798 


0.168±0.004 


214.7±2.0 


6 


50499 


0.193±0.002 


206.5±0.2 


7 


50810 


0.181±0.004 


205.9±0.8 


7 


51017 


0.142±0.005 


221.0±2.0 


7 


51017 


0.147±0.006 


223.0±2.0 


7 



Notes: Observation dates are listed as offsets relative to 
Julian Day 2400000. 

References: 1) Mariotti et al., 1990; 2) Henry & Mc Carthy 
1993; 3) This paper, Kitt Peak, 2D speckle imager; 4) This 
paper, ESO 3.6m, COME-ON adaptive optics system; 5) 
CFHT, CIRCUS 2D speckle imager; 6) ESO 3.6m, ADONIS 
adaptive optics system; 7) This paper, CFHT, PUE'O 
adaptive optics system. 

ror commands, as described by Veran et al. (1997). For 
pre-1997 measurements it was instead obtained from ob- 
servations of a reference single star of similar R-band mag- 
nitude. Astrometric calibration fields such as the central 
region of the Trapezium Cluster in the Orion Nebula (Mc- 
Caughrean & Stauffer 1994), were observed to accurately 
determine the actual detector plate scale and position an- 
gle (P.A.) origin. Uncertainties on these parameters do 
not apreciably contribute to the overall separation error 
for the small separations in the Gl 570BC system. 

The separation, position angle and magnitude differ- 
ence between the two stars were determined using uv plane 
model fitting in the G1LDAS (Grenoble Image and Line 
Data Analysis System) software, as well as with the de- 
convolution algorithm described by Veran et al. (1999), 
coded within IDL. With approximate initial values of the 
positions of the two components along with the PSF ref- 
erence image, the fitting procedures gave as output the 
flux and pixel coordinates of the primary and secondary. 
Application of the astrometric calibrations then yields the 
desired parameters. 



2.2. Radial velocity measurements 

All radial velocity measurements are listed in Table 2 
(available in electronic form only). Most of them were ob- 
tained with the two CORAVEL radial velocity scanners 
(Baranne et al. 1979) on the Swiss lm telescope at Ob- 
servatoire de Haute Provence (France) and on the Dan- 
ish 1.54m telescope at La Silla (Chile). The earlier data 
were previously published by Duquennoy & Mayor (1988) 
and the system has been regularly observed since then. 
Gl 570BC is relatively bright for the CORAVEL instru- 
ments (V=8.09, Leggett 1992), but both stars are poor 
matches to the fixed K0III correlation mask. They there- 
fore only produce shallow correlation dips (of which exam- 
ples were displayed in Duquennoy & Mayor (1988)) and 
their velocities are measured with typical precisions of re- 
spectively 0.7 and 4 km/s, instead of the usual CORAVEL 
precision of 0.3 km/s. The measurements of the M3V sec- 
ondary in particular were at the limit of the CORAVEL 
capabilities. During the orbit adjustment discussed be- 
low they were found to have sizeable systematic errors 
at phases where the two profiles are even slightly blended. 
This didn't measurably affect the derived orbital elements 
since these noisier measurements carried essentially no 
weight anyway, but they were nonetheless ignored in the 
final solution. For clarity, they are also not plotted in Fig- 
ure 1. 

Over the last three years, we have obtained consider- 
ably more accurate measurements with the ELODIE spec- 
trograph (Baranne et al. 1996) on the 1.93m telescope 
of Observatoire de Haute Provence and the CORALIE 
(Queloz et al. in preparation) spectrograph on the Swiss 
1.2m Euler telescope at La Silla (Chile). These echelle 
spectra were analysed by numerical cross-correlation with 
an M4V one-bit (i.e. 0/1) mask, as described by Delfosse 
et al. (1999b). Radial velocities were initially determined 
by adjusting double gaussians to the correlation profiles. 
Those however had systematic phase-dependent residu- 
als during the orbital adjustment, which were particularly 
large (400 m/s rms) for the fainter secundary star. Upon 
analysis, they were found to originate from the low level 
wiggles in the correlation profile of the primary: while the 
core of this profile is well described by a gaussian func- 
tion, its baseline doesn't drop to zero as a gaussian would, 
but instead keeps oscillating at the ±0.2% level, about 
one tenth of the depth (2%) of the secondary star's dip. 
Depending on its position on this oscillating background, 
the measured velocity of the fainter star could therefore 
be incorrect by up to 1 km/s. Similar errors of course also 
affected the measured velocities of the primary star, but 
they were smaller by the square of the relative depths of 
the two correlation dips, or a factor of 16. Thanks in part 
to the excellent stability of the spectrograph, the shape 
of the correlation profile for a given star is very stable. 
We could therefore obtain an excellent estimate of the 
wings of the intrinsic correlation profile of each star, by 
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averaging all profiles of the system, after aligning them at 
the measured velocity of the star and blanking all pixels 
within two profile widths of the velocity of the other star. 
The residuals of a gaussian adjustment to these average 
profiles are then subtracted from all correlation profiles, at 
the measured velocity of each star. This decreases the fluc- 
tuation level in the profile baseline to a level of ±0.05%. 
The radial velocities measured by a double gaussian fit to 
these corrected profiles have typical accuracies of 40 m/s 
for the primary and 100 m/s for the fainter secondary. 
These residuals are twice larger than would be measured 
for single- lined spectra with equivalent S /N ratios and cor- 
relation dip parameters, but show no systematic phase de- 
pendency. They should thus cause no systematic errors on 
measured parameters. 

3. Orbit adjustment 

3.1. The adjustment code 

The program used for the orbit adjustment, ORBIT, de- 
rives from the code of Tokovinin (1992). Like its progeni- 
tor, it performs a least square adjustment to all available 
spectroscopic and "visual" observations, with weights in- 
versely proportional to the square of their standard errors. 
Besides a number of cosmetic changes, the modifications 
and additions to Tokovinin (1992) include allowance for 
multiple radial velocity zero-points, the use of the more 
robust Levenberg-Marquardt minimisation algorithm, as 
well as direct support for triple systems. Standard errors 
for derived parameters (masses, parallax, etc) are com- 
puted from the full covariance matrix of the linearized 
least square adjustment, rather than from just the stan- 
dard errors of the orbital elements (i.e. the diagonal terms 
of the covariance matrix). In addition ORBIT can es- 
timate confidence intervals for both derived parameters 
and orbital elements through Monte-Carlo experiments. 
For well constrained binaries (like Gl 570BC) these inter- 
vals are consistent with the analytic standard errors and 
gaussian statistics, incidentally providing a check against 
gross errors in the analytic gradients. For noisier systems 
though, the vicinity of the minimum of the x 2 surface be- 
comes significantly non-quadratic. The error statistics for 
the parameters are then substantially non-gaussian and 
asymetric, and the Monte-Carlo confidence intervals be- 
come an essential feature. 

ORBIT also accepts some additional data types, be- 
yond spectroscopic velocities and "visual" (p, 9) or (X,Y) 
pairs. These additions, all of which are used in the present 
paper, include: 

— ID projected separation measurements, produced by 
early IR scanning speckle observations and lunar oc- 
cultations; 

— parallaxes, which can be seen as an additional link 
between the velocity amplitudes, the eccentricity, the 
period and the semi-major axis; they enter in the 



least square adjustment with their standard error, pro- 
viding an optimal unified description between orbital 
parallaxes of SB2+visual pairs and mass ratios from 
SBl+visual pairs with a known parallax; 
— cross-correlation "dips", produced by correlation ve- 
locimeters and cross-correlation analysis of echelle 
spectra; to support these profiles ORBIT complements 
the usual list of 10 orbital elements by the width and 
the equivalent width of one gaussian profile for each 
spectroscopic component, and adjusts the 10 orbital 
elements and the profile parameters directly to the en- 
semble of all cross-correlations, rho-theta pairs, par- 
allaxes, etc. In a similar manner, modern analyses of 
long baseline interferomctric data sometimes adjust or- 
bital elements (together with magnitude differences) 
directly to the uv data, avoiding the intermediate step 
of p, 8 extraction (Hummel & Amstrong 1992), as do 
adjustments of visual orbits to Hipparcos Transit Data 
(Quist & Lindegren 1999; Soderhjelm 1999). A practi- 
cal inconvenient of this bypassing of intermediate steps 
is that no individual radial velocities (or p, 8, or separa- 
tions on great circles) are available for compact pub- 
lication, and that the usual figures of radial velocity 
curve (or visual orbit) loose sense and become illus- 
trative only. On the positive side however, direct or- 
bital adjustment to the profiles greatly decreases the 
effective number of degrees of freedom, compared with 
orbital adjustment to velocities measured from indi- 
vidual profiles. We found that this often improves the 
error bars on the spectroscopic orbital elements by a 
factor of two, or better. Global adjustment also consid- 
erably reduces the susceptibility of the orbital solution 
to "pulling" of the velocity of the weaker component 
towards and away from that of its brighter compan- 
ion for profiles which are only imperfectly gaussian. In 
addition, and to further diminish this sensitivity, we 
produce an average recentered profile for each com- 
ponent, determine residuals to the best fitting gaus- 
sian, and subtract those residuals from the individual 
measurements to then determine iteratively improved 
parameters. 

ORBIT runs on Unix systems and can be obtained from 
the authors upon request. 

To illustrate our methods and the power of combining 
different data types for the same binary, as probably first 
advocated by Morbey (1975), we present 3 distinct orbits 
for Gl 570BC, incorporating successively more of the in- 
formation available to us. The speckle and adaptive optics 
data alone turn out to be insufficient to properly define 
a visual orbit, largely because they never resolve this ec- 
centric system close to its periastron (Fig. 2), where the 
minimal separation is only ^12 milliarcsecond. As a con- 
sequence these data alone leave indeterminate a combina- 
tion of inclination and semi-major axis with eccentricity. 
The most restrictive orbit we can present in Table |^ is 
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therefore a spectroscopic one. We then present a combined 
spectroscopic and visual orbit, and finally an orbit which 
in addition uses the independent trigonometric parallax of 
the Gl 570 system. 

The longitude of the periastron is given with the 
spectroscopic convention, and thus refers to the primary. 
180°must be added to uj to obtain the visual convention. 
Table |l| lists the individual speckle and adaptive optics 
measurements and Table 2 (only available electronically) 
gives the same information for the radial velocities. Note, 
though, that the final orbit was directly adjusted to the 
individual cross-correlation profiles, while the velocities 
given in Table 2 were determined separately by adjust- 
ing two Gaussian curves to each profile. These velocities 
are given here mostly for illustrative purposes, as plotted 
in Fig. 1. Their use would give a slightly different (and 
noisier) solution, but may nonetheless prove convenient. 
There is no compact way to publish the full information 
needed to reproduce our analysis, since it does not proceed 
through determining radial velocities as an intermediate 
step. We will make the digital correlation profiles available 
upon request to the first author. Fig. 2 shows the visual 
data and orbit. 

3.2. Spectroscopic orbit 

As mentioned above, the spectroscopic orbit presented 
in Table || was adjusted directly to the correlation pro- 
files rather than to the extracted radial velocities. As can 
be seen, the spectroscopic quantities Mxsin 3 i are deter- 
mined with very high accuracies of 0.2%. The relatively 
few ELODIE measurements contribute considerably to 
the overall precision of the orbital solution, and ignoring 
them would for instance degrade the standard errors of 
the Mxsin 3 i by an order of magnitude. The much more 
numerous CORAVEL measurements of the primary star 
by contrast only contribute to an improved orbital pe- 
riod, thanks to their much longer timespan. They other- 
wise carry very little weight in the solution. 

3.3. Spectroscopic+visual orbit 

Inclusion of the speckle and adaptive optics data in the 
adjustment leaves all spectroscopic elements essentially 
unchanged (Table ||), but determines the three otherwise 
unknown orbital elements: the semi-major axis (a), the 
inclination (i), and the orientation of the projected or- 
bit on the sky (CI). Individual stellar masses as well as 
an orbital parallax for the system can then be derived 
from the full complement of orbital elements, and are also 
listed in Table ||. The standard errors computed from the 
covariance matrix were checked through Monte-Carlo sim- 
ulations, which resulted in confidence intervals that were 
fully consistent with gaussian errors of the stated disper- 
sion. This verifies that non-linearities in the least square 
adjustment are negligible for this very well constrained 
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Fig. 1. Radial velocity orbit of the Gl 570BC system. 
Filled triangles represent measurements of the primary 
star, and filled squares measurements of the secundary 
stars. The ELODIE mesurements can be distinguished 
from the CORAVEL ones by their much smaller error 
bars. 



system. The two masses are determined here with 1.6% 
accuracy, and the orbital parallax to within 2.2 milliar- 
second (1.3%). 

The availability of Hipparcos parallaxes represents an 
opportunity to independently verify the orbital parallax, 
and thus to globally check the orbital solution for system- 
atic errors. Though somewhat noisier than typical for a 
V=5.7 star, the Hipparcos catalog parallax of Gl 570A 
is well determined, tt=0.1693±0.0018". The parallax for 
Gl 570BC itself has very large error bars (ct(tt)=0.033") 
in the Hipparcos catalog (ESA 1997), because the un- 
accounted orbital motion with P^0.8 year strongly cou- 
ples into the parallax solution over the limited lifetime of 
the Hipparcos satellite. Fortunately, Soderhjelm (1999) re- 
cently reanalysed the Hipparcos intermediate transit data, 
accounting for the orbital motion within Gl 570BC, and 
obtained sharply reduced error bars for the trigonometric 
parallax: tt=0.1697"±0.0010". The two astrometric deter- 
minations are mutually consistent and agree with the or- 
bital parallax of 0.1710±0.0022", to within better than 
1 a. 

3.4. Spectroscopic+Visual+Parallax orbit 

Instead of using the independent trigonometric parallaxes 
as a sanity check for the orbital solution, ORBIT offers the 
option to consider it as an additional observation, linking 
i, e, P, a and K1+K2. It is then included in the combined 
least square adjustment, together with the radial velocities 
and angular separations. This very significantly improves 
the determinacy of the least square system, and in par- 
ticular reduces the standard errors of the semi-major axis 
by a factor of 2. We note here that it would be prefer- 
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Table 3. Orbital elements and derived parameters 



Element 

V (km/s) 
P (days) 
To (JD) 
e 

a (arcsec) 
fii (deg) 
Ui (deg) 
i (deg) 
Ki (km/s) 
K 2 (km/s) 
Wi (km/s) 
W 2 (km/s) 
EWi (km/s) 
EW 2 (km/s) 

Derived parameters: 
Mi xsm 3 i(M ) 
M 2 xsm 3 j(M ) 
Mi (M ) 
M 2 (M ) 
7t (arcsec) 



Spectro 

28.664±0.008 
308.884±0.004 
270.220±0.011 
0.7559±0.0002 



127.56±0.05 

18.187±0.008 
27.325±0.026 

6.169±0.008 

5.643±0.030 
0.4532±0.0005 
0.1041±0.0005 

-0.633±0.065 



0.5082±0.0011 
0.3383±0.0005 



Spectro+OA Spectro+OA+7r 



28.665±0.008 
308.884±0.004 
270.220±0.011 
0.7559±0.0002 
0.1514±0.0015 
196.2±0.7 

127.56±0.05 
107.2±1.0 

18.187±0.008 
27.324±0.026 

6.169±0.008 

5.642±0.030 
0.4532±0.0005 
0.1041±0.0005 

-0.633±0.065 



0.5082±0.0011 
0.3383±0.0005 
0.583±0.009 
0.388±0.006 
0.1710±0.0022 



28.665±0.008 
308.884±0.004 
270.220±0.0110 
0.7559±0.0002 
0.1507±0.0007 
195.9±0.5 

127.56±0.05 
107.6±0.7 

18.187±0.008 
27.325±0.026 

6.169±0.008 

5.642±0.030 
0.4532±0.0005 
0.l04l±0.0Cg5 1 

-0.633±0.065 



0.5083±0.0011 
0.3383±0.0005 
0.586±0.007 
0.390±0.005 
0.1698±0.0009 



Notes: Epochs are listed as offsets relative to Julian Day 
2450000. Wi and W 2 are the full widths to half power of the 
Gaussian function for the two components, while EWi and 
and EW 2 are their equivalent widths. The other orbital 
elements have their usual meaning. 








-0.05 



-0.1 



Fig. 2. Visual orbit of the Gl 570BC system. Separations 
are in arc seconds, North is up and East is left. 



able, at least conceptually, to adjust an orbit directly to 
the Hipparcos transit data rather than to the Hipparcos 
parallax. ORBIT does not yet support this data type how- 
ever, which we plan to add in the near future. 

The contributions to the overall \ 2 of the different data 
types included in the solution are approximately consis- 
tent with their respective number of measurements. This 
indicates that there are no large systematic errors in any 
one data type, and that the adopted standard errors are 
at least approximately correct. This orbit (last column of 
Table ||) has the smallest errorbars (and smallest covari- 
ances) for all orbital elements. It is consistent with all pre- 
viously published orbits to within their stated error bars, 
after allowing for the 180°ambiguity in the identification 
of the ascending node (i.e. Q) in purely visual orbits. We 
adopt it as our preferred solution for the rest of the discus- 
sion. The inclusion of the trigonometric parallax informa- 
tion improves the relative accuracy of the two masses to 
1.2%: Mi = O.586±O.OO7M and M 2 = 0.390±0.005 M Q . 
They are consistent with the values obtained by Mariotti 
et al. (1990) to within a fraction of their quoted stan- 
dard deviation, as well as with the mass sum derived by 
Soderhjelm (1999) from an analysis of the Hipparcos tran- 
sit data, but improve on their accuracy by over a factor 
of 5. These accuracies are among the best obtained to 



date for non-eclipsing binaries, but still do not match the 
~0.2% obtained for the spectroscopic Mxsin 3 i. There is 
thus still room for improvement in the astrometric mea- 
surements. We are conducting some long baseline inter- 
ferometric observations of Gl 570BC (Segransan et al. in 
progress) with IOTA (Millan-Gabet et al. 1999), which 
have the potential to ultimately match the sub-% spectro- 
scopic accuracy for the masses. The "resulting" parallax in 
the last column of Table |), 0.1698±0.0009", is computed 
as an orbital parallax from the elements of this combined 
orbit, which themselves take the trigonometric parallax 
into account. It is an optimal combination of the distance 
information available in the trigonometric parallax and in 
the orbit. In the present case is mostly determined by the 
weighted average of the two HIPPARCOS parallaxes, with 
little contribution from the orbit. 

4. Luminosities and colours 

Photometric measurements of the combined light of 
the system have been made by a number of au- 
thors over a broad range of colours. They are sum- 
marised and homogenized in the critical compilation of 
Leggett (1992), from which we adopt: U=10.79, B=9.57 
,V=8.09, R=7.09, 1=5.97, J=4.75, H=4.14, K=3.93, 
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Table 4. Absolute magnitudes of the two components 





Primary 


Secondary 


Mj 
M„ 
Mk 
M l 


6.213±0.033 
5.613±0.033 
5.396±0.033 
5.167±0.111 


7.403 ±0.039 
6.763 ±0.039 
6.576 ±0.039 
6.227 ±0.159 



L=3.77, L'=3.67, with quoted errors of <3% on IJHK. The 
system has colours consistent with those of solar metallic- 
ity stars (Lcggett 1992), as expected from the spectro- 
scopic metallicity discussed below. 

The difference in brightness between the components 
of Gl 570BC has been measured on several occasions 
in different bandpasses. Mariotti et al. (1990) obtained 
magnitude differences of A, 7 =1.27±0.12, A JJ =1.19±0.12, 
A^=1.12±0.07, and A' L =1.06±0.17 from their ID scan- 
ning speckle observations. HMcC published another esti- 
mate in the J band, Aj=1.30±0.04. Our own adaptive 
optics measurements provide Ak =1.18±0.04, as well as 
A m =1.18±0.04 in a narrow band Br 7 (2.166/xra) filter. 
We estimate a correction of Ak — As r =0.04, using syn- 
thetic photometry generated from NextGen (Hauschild, 
Allard, & Baron 1999) model spectra of effective temper- 
atures that bracket those apropriate for M dwarfs of the 
luminosities of the two components of Gl 570BC. This 
provides a second determination of the K band flux ratio, 
Atf=1.22±0.04. 

While those measurements are mutually consistent for 
every filter, their run with wavelenghth is only marginally 
compatible with the known IR colours of M dwarfs 
(Leggett 1992). The J—H and (particularly) J—K colours 
of early and mid-M dwarfs define a remarkably flat plateau 
(Lcggett 1992): when the spectral type of solar metallic- 
ity stars (apropriate for Gl 570BC) varies between MOV 
(M K -4.8) and M5.5V (M K -7.9), J - H only changes 
from 0.68 to 0.57, and J-K just from 0.85 to 0.87 (Leggett 
1992) The two components of Gl 570BC respectively have 
M K -5.4 and M K -6.6 and their J - H and J-K 
colours are therefore expected to only differ by about - 
0.04 and ±0.01. Aj should therefore be almost identical 
to A if, whereas the observations give Aj-Ak =0.14±0.04, 
and Aj-Aii should be -0.04, whereas the observations in- 
dicate Aj-Aif=±0.11±0.12. As the integrated colours of 
the system match the expected values, this inconsistency 
must come from some of the measured magnitude differ- 
ences. We suspect that it traces back to an overestimated 
contrast in the J and H band observations of both Mari- 
otti et al. (1990) and HMcC, since speckle techniques have 
a known bias towards underestimating the relative flux of 
faint components (Perrier 1988). The K band adaptive op- 
tics observations are in principle immune to this bias, and 
the K band speckle observations of Mariotti et al. (1990) 
should be relatively unaffected, thanks to the low D/r 
(where D is the telescope diameter and tq is the Fried 
parameter of the atmosphere) at this longer wavelength. 



We therefore tentatively adopt as the basis of our mag- 
nitude difference system the mean of the three K band 
measurements, Ak =1.18±0.03. From this value we then 
derive preferred values of Aj=1.19 and A/i=1.15, but we 
will also consider the published J and K flux ratios as an 
alternative. This discrepancy contributes significant un- 
certainty to the analysis, and better measurements of the 
flux ratios at J and H would be of considerable interest. 
To date all our adaptive optics measurements in the J 
and H bands have unfortunately been obtained at phases 
when the secondary star overlaps the first Airy ring of the 
primary for these wavelengths. These circumstances maxi- 
mize the uncertainties in differential photometry from par- 
tially corrected adaptive optics images (Veran et al. 1999), 
and these data therefore contribute no useful information 
on the flux density ratio. 

Absolute magnitudes were derived from the parallax 
and the apparent magnitudes of the individual stars. The 
absolute magnitudes of the brighter star have uncertain- 
ties which are dominated by those of the photometry of the 
system, while those for the secondary have some contribu- 
tion from the magnitude difference. The parallax doesn't 
appreciably contribute in either case. 

5. Discussion and conclusions 

The accurate masses and absolute magnitudes that we 
have obtained for the Gl 570BC system represent a new 
benchmark for model calculations (e.g. Baraffe et al. 
1998) and an independent check of the empirical mass- 
luminosity relations (HMcC). The constraints which they 
bring to the models are largely complementary to those 
coming from the eclipsing binaries, whose absolute radii 
can be determined very accurately but whose larger dis- 
tances on the other hand contribute significant uncertain- 
ties to the absolute magnitudes. 

As long emphasized by theoreticians, and by observers 
of more massive stars (e.g. Andersen 1991), there is 
however no such thing as one single mass-luminosity 
relation: stellar luminosities depend on chemical com- 
position as well as on mass (in general they depend on 
age too, though not in the age and mass range discussed 
here). Quantitative metallicity determinations however 
are notoriously difficult for M dwarfs (e.g. Viti et al. 
1997; Valenti et al. 1998). Observers in this field usually 
have to resort to photometric metallicity estimators 
(Lcggett 1992) which are only approximately calibrated, 
or otherwise assume by default a solar metallicity. Thanks 
to its physical association with the hotter Gl 570 A (K4V), 
the Gl 570BC pair represents a rare case of two M dwarfs 
with a very well determined spectroscopic metallicity. 
Its accurate masses are thus fortunately matched with 
excellent metallicities. Hearnshaw (1976) first measured 
the metallicity of Gl 570A from high resolution elec- 
tronographic spectra and obtained [Fe/H]=+0.01±0.15. 
More recently, Feltzing & Gustafsson (1998) measured 
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Fig. 3. Near-IR Mass-Luminosity diagrams for mid-M 
dwarfs. The square symbols with thick error bars corre- 
spond to the two components of Gl 570BC, and the thin 
error-bars represent data from HMcC. The solid and dot- 
dashed lines respectively represent the empirical HMcC 
analytical Mass-Luminosity relation and the theoretical 
[M/H]=0.0 10 Gyr isochrone of Baraffe et al. (1998). 

[Fe/H]=0.04±0.02(random)±0.1(systematic) from high 
S/N 

R=10 5 echelle CCD spectra. These authors find some 
evidence for NLTE overionisation of Fe into Fe + , but the 
derived elemental Fe abundance is unaffected, as Fe is 
overwhelmingly neutral in the photosphere of a K4 dwarf. 
Quite conveniently for comparison with published models, 
the Gl 570 system thus has a truly solar metallicity. 

Figure ^ compares the mass and luminosity of the two 
components of Gl 570BC with the Baraffe et al. (1998) 
10 Gyr solar-metallicity isochrone. These models consis- 
tently combine stellar evolution models (e.g. Chabrier & 
Baraffe 1997) and non-grey atmospheric models (Allard & 
Hauschild 1995; Hauschild ct al. 1999). The present gen- 
eration of these evolutionary models still uses non-dusty 
atmospheres, but dust only becomes relevant at effective 
temperatures significantly lower than those of Gl 570BC 
(Allard 1998). The Baraffe et al. models are consistently 
slightly brighter in all 3 bands than the two stars, by 0.08 
to 0.15 magnitude. While this level of agreement is al- 



ready very comforting, the discrepancies are significant 
at the ~3cr levels and may point towards remaining low 
level deficiencies of the theoretical models. If one adopts 
the measured J and H band flux ratios, rather than the 
smaller values that we have estimated from the K band 
measurement, the agreement with the models improves 
only slightly for the primary star, as its magnitude only 
weakly depends on the exact value of the large flux ratio. 
At the same time this choice degrades the agreement for 
the fainter star by a larger factor, and the overall agree- 
ment with the models is significantly worse. 

Figure || also shows the data points of HMcC, as well 
as their analytic representation of those data. The agree- 
ment is essentially perfect with the J band HMcC relation, 
while the H and K band relations are slightly discrepant, 
by respectively 0.1 and 0.15 magnitudes. We note that 
the HMcC mass-luminosity relations are only consistent 
with the empirical M dwarf colours of Leggett (1992) at 
this 0.10-0.15 magnitude level, even though the HMcC 
photometry mostly traces back to Leggett (1992). This 
is because the HMcC mass-magnitude relations for J, H 
and K were adjusted independently, without explicit forc- 
ing of colour consistency. The perfect agreement with the 
J band HMcC fit is therefore probably fortuitous to some 
extent, and the 0.10 to 0.15 discrepancy for the H and K 
bands probably represents a more realistic estimate of the 
accuracy of those analytic fits around 0.5 M©. 
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± 


.144 


50839.715 


37.627 ± .033 


14.939 


± 


.140 


50851.704 


38.237 ± .031 


14.096 


± 


.117 


50852.707 


38.283 ± .035 


14.000 


± 


.126 


50853.704 


38.273 ± .033 


13.919 


± 


.127 


50854.704 


38.309 ± .030 


13.906 


± 


.121 


50857.713 


38.417 ± .033 


13.769 


± 


.129 


50886.622 


15.867 ± .023 


47.349 


± 


.086 


50889.637 


5.941 ± .036 


62.600 


± 


.125 


50890.645 


3.907 ± .033 


65.285 


± 


.131 



Julian Day 




v 2 




2400000+ 


(in km/s) 


(in km/s) 




CORALIE measurement 


51197.874 


7.546 ± .030 


60.549 ± 


120 


51198.870 


5.057 ± .030 


63.986 ± 


130 


51200.880 


2.391 ± .030 


68.186 ± 


120 


51207.859 


5.868 ± .030 


62.853 ± 


120 


51209.886 


7.488 ± .030 


60.661 ± 


130 


51211.778 


8.883 ± .030 


58.694 ± 


120 



